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OPERATING  CHARACTERISTICS  FOR  DETECTION  OF  A 
FADING  SIGNAL  IN  M  ALTERNATIVE  LOCATIONS 
WITH  D-FOLD  DIVERSITY 

INTRODUCTION 

Before  the  transfer  of  information  between  a  transmitter  and  receiver  can 
take  place,  synchronization  is  necessary;  that  is,  the  receiver  must  ascertain 
the  time  delay  and  frequency  shift  of  the  received  signal  before  correct  deci¬ 
sions  about  the  information  content  of  an  ensuing  message  can  be  made.  For 
large  a  priori  uncertainties  about  the  transmitter's  range  and  relative  velocity, 
the  receiver  must  conduct  a  search  —  during  this  alert  phase  —  of  the  numerous 
possible  locations  (time-delay  frequency-shift  cells)  to  determine  the  precise 
time  delay  and  frequency  shift  of  the  received  signal,  and  to  determine,  in  fact, 
whether  a  signal  is  present  at  all.  Furthermore,  in  order  to  combat  the  possi¬ 
bility  of  signal  fading,  the  transmitter  may  employ  diversity  by  dividing  the 
available  signal  energy  into  a  number  of  branches.  The  receiver  must  then  know 
the  order  and  pattern  of  diversity  and  combine  the  energies  of  the  appropriate 
diversity  branches  before  making  a  decision  of  absence  or  presence  and  loca¬ 
tion  of  a  possible  signal.  The  number  of  locations  that  the  received  signal  may 
occupy,  if  present,  is  denoted  by  M;  the  order  of  diversity  employed  is  denoted 
by  D.  The  objective  of  this  report  is  to  determine  the  performance  (in  terms 
of  appropriate  probabilities  of  various  types  of  errors)  of  the  synchronization 
procedure  as  a  function  of  M,  D  and  the  received  signal  and  noise  power  levels. 
Additionally,  the  optimum  order  of  diversity  and  optimum  fractionalization  of 
available  signal  energy  for  specified  performance  are  to  be  determined. 

If  the  arrival  angle  of  the  potential  signal  is  not  known  to  the  receiver,  a 
search  must  be  conducted  in  this  variable  as  well  as  in  time  delay  and  frequency 
shift.  If  the  parameter  M  is  increased  to  include  this  additional  uncertainty, 
the  case  is  subsumed  under  the  earlier  framework,  where  M  denotes  the  num¬ 
ber  of  potential  signal  locations,  whether  they  be  cells  in  time  delay,  frequency 
shift,  or  arrival  angle.  More  generally,  M  is  simply  the  number  of  locations 
in  which  a  signal  can  be  found,,  if  present,  regardless  of  the  physical  cause  of 
uncertainty. 

Although  the  discussion  has  been  couched  in  a  detection  context,  the  frame¬ 
work  also  covers  multiple  alternative  communication  with  thresholding.  For 
example,  suppose  that  either  no  signal  is  transmitted  or  one  of  M  signals  is 
transmitted.  Then,  for  fading  signals  and  D-fold  diversity,  the  results  to  be 
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presented  are  directly  applicable  to  determining  the  error  probabilities.  As 
particular  cases,  M  -  1  corresponds  to  on-off  (binary)  communication  with 
diversity;  M  -  2  and  a  zero  threshold  corresponds  to  binary  frequency-shift¬ 
keying  with  diversity;  and  M  -•  2  and  a  zero  threshold  corresponds  to  M-ary 
communication  with  diversity,  as  previously  analyzed1  and  numerically  evalu¬ 
ated.  2  For  M  =  1,  some  results  on  tills  problem  are  available3  (when  D  here 
is  identified  with  M  there). 

To  develop  the  subject,  two  types  of  probabilities  of  mistakes  are  defined, 
and  generic  expressions  for  their  values  are  presented.  Then  three  different 
signal  and  noise  models  —  all  of  which  yield  the  same  statistics  for  the  observed 
variables  —  are  described,  and  the  optimum  processor  is  derived.  Finally,  the 
actual  evaluation  of  the  correct-decision  and  false-alarm  probabilities  is  under¬ 
taken  and  plotted,  and  relevant  observations  and  conclusions  are  extracted  from 
the  numerical  results. 

DEFINITIONS  OF  ERROR  PROBABILITIES 

For  a  signal  that  can  be  found  in  one  of  M  locations,  let  xp  •  •  •  ■  be 
the  M  decision  variables  upon  which  a  decision  must  be  reached  as  to  whether 
signal  is  absent  or  whether  a  particular  signal  is  present.  If  D-fold  diversity 
is  employed,  the  vanable  xk  is  an  appropriate  combination  (as  yet  undefined) 
of  the  D  branch  outputs  utilized  by  signal  number  k.  If  threshold  A  is  uti¬ 
lized,  it  is  decided  that 

signal  is  absent  if  max  (x-, , . . . ,  xM)  <  A 

(1) 

signal  k  is  present  if  max  (Xp  . . . ,  xM)  =  x^  >  A 

Then  the  probability  of  false-alarm  is  given  by 

PFA  =  1  -  Prob  jmax  (xp ...»  XM)  <  A  |N  |  ,  (2) 

where  |N  denotes  conditioning  on  noise-alone  present  (that  is,  no  signal  pres¬ 
ent).  In  words,  (2)  is  the  probability  that  one  or  more  of  the  M  decision  vari¬ 
ables  exceeds  the  threshold  for  noise-alone  present. 

When  a  signal  is  present,  two  types  of  mistakes  can  occur.  In  order  to 
express  these  quantities,  without  loss  of  generality,  let  signal  number  1  be 
transmitted.  The  probability  of  a  missed  decision  is  defined  as 

pMD  ~  Prob  imax  <xl»  *  •  *  *  XM>  <  A  ls+  N  1  > 
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where  |S*  N  denotes  conditioning  on  signal-plus-noise  present.  That  is,  it  is 
erroneously  decided  that  no  signal  is  present  because  the  threshold  was  not  ex¬ 
ceeded. 

On  the  other  hand,  the  probability  of  an  incorrect  decision  is  defined  as 
M 

PID=  ^  Prob  jmax  (Xi, . ,  xM)  -  xk  >A|S+N  |  .  (4) 

k-2 

That  is,  one  of  the  wrong  decision  variables  dominates  and  exceeds  the  thresh¬ 
old,  thereby  leading  to  the  wrong  decision  about  which  signal  was  transmitted.  * 

Finally,  the  probability  of  a  correct  decision  is  defined  as 

PCD  =  Prob  {max  (x^  . . . ,  xM)  =  x1  >A  |S+  N}  .  (5) 

This  is  the  probability  that  the  correct  decision  variable  is  largest  and  exceeds 
the  threshold.  Since  one  of  the  events  described  in  (3),  (4),  or  (5)  must  happen 
for  signal  present,  we  have 


PMD  + 


P  + 
ID 


PCD  1 


An  alternative  form  of  (G)  is 


(6) 


P 


MD  + 


(7) 


where  the  left  side  of  (7)  is  the  sum  of  the  probabilities  of  the  two  types  of  mis¬ 
takes  when  signal  is  present. 


An  average  error  probability  could  be  defined  (if  desired)  in  the  following 
manner:  an  error,  for  a  signal  transmitted,  occurs  with  probability 


Pe  <S  ,  N>=PMDtPID'’1-pCD  ' 


An  error,  for  no  signal  transmitted,  occurs  with  probability 


P 

e 


(N)  =  P 


FA  ’ 


<8) 

(9) 


*An  on-off  communications  system  corresponds  to  M  =  1;  in  this  case, 
Pjq  is  not  defined  (that  is,  it  has  no  meaning). 
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Then  the  average  error  probability  is 

Pe  -  Pe(S+  N)  Prob  |S  +  N|  +  Pe(N)  Prob  |n| 

(10) 

=  (l~pCD)Prob  iS+Nl  ■»  p FA  Prob  | N |  , 
where  Prob  jS+Nj  and  Prob  |n|  are  the  a  prion  probabilities. 
INTEGRAL  EXPRESSIONS  FOR  PROBABILITIES 


Analytic  evaluation  of  the  probabilities  in  (2)  through  (5)  is  generally  im¬ 
possible  unless  the  variables  xA? . . . ,  xpj  arc  all  statistically  independent  of 
each  other.  From  this  point  on,  this  crucial  assumption  of  independence  will  be 
adopted.  It  is  also  assumed  that  the  random  variables  jxj.}  not  containing  a 
signal  have  identical  probability  density'  functions  pQ(  . ).  The  random  variable 
containing  a  signal  will  have  probability  density  function  pj{  • ).  The  cumulative 
distribution  functions  are  denoted  by  P0(- )  and  Pi(*),  respectively. 

The  quantities  in  (2)  through  (5)  then  follow  immediately: 


PFA  1 


r  rA  iM  r  iM 

-  I  dx  po(x)  =  1  -  p0(A) 


rA  r 

rA  I*1"1 

PMD=  /  ^  P1(X) 

/  dy  p0(y)  =  pj(A) 

J-*  L 

J-oc  J 

ID 


['*o«A)l 

»  (M-  I)  dx  pQ(x)  J  dy  Pj(y)  f  j  dz  p0(z)l 

r  ~]M  -  2 

=  (M-l )J  dx  pQ(x)  Pi(*)  |p0(J 

A 


M-2 


(ID 

(12) 


(lo) 


and 


r 

P _ -  J  dx  pA(x) 


CD 


[£dy  P0«] 


M-l  /— > 

-  /  dx 

Jk 


PA(x) 


po<x> 


M-l 


.  (14) 
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(As  a  check,  notice  that 
<1 


f  <1  I  r  tm"1i  T  1M'1 

%  ‘  To  i*  * jV* ['*.«  | =  '  -  'YA> L'o<a>J  ■ 


-1-P 


MD 


<15> 


which  agrees  with  (0). ) 


Of  the  four  proh;ibilities  defined  in  (2)  through  (5),  onlythree  are  independ¬ 
ent;  see  ((J).  Furthermore,  Pp^  imd  PMD  can  be  found  very  from  l11) 

and  (12)  once  the  cumulative  probability  distributions  P0(- )  and  Pj(- )  are 
known.  If  (M)  is  numerically  integrated  in  order  to  determine  P^p),  PlD  is 
immediately  available  by  use  of  (6).  The  discussion,  then,  will  be  concentrated 
on  the  two  probabilities  PpA  and  PCD>  these  are  also  the  only  two  probabilities 
needed  to  evaluate  the  average  error  probability  defined  in  (10).  Therefore, 
only  one  numerical  integration,  that  in  (14),  need  be  conducted.  An  error 
anal}  sis  for  (14)  and  some  interrelationships  of  the  probabilities  in  (11)  through 
(14)  are  presented  in  appendix  A. 


SIGNA  I.  AND  NOISE  MODELS  AND  THE  OPTIMUM  PROCESSOR 


In  the  previous  section,  the  M  decision  variables  x^, ...»  xm  were  pre¬ 
sumed  to  be  available  for  comparison  with  a  threshold.  Now  w'e  will  backtrack 
to  ascertain  how  these  decision  variables  can  be  determined  in  the  first  place. 
The  first  model  to  consider  is  that  in  which  the  received  signal  consists  of  nar¬ 
rowband  Gaussian  processes  on  D  statistically  independent  branches.  The 
particular  D  branches  occupied  by  the  transmission  of  one  of  the  M  alterna¬ 
tives  do  not  overlap  any  of  the  (M  -  1)  D  branches  utilized  for  the  other  signal 
alternatives.  All  MD  branches  are  subject  to  additive  narrowband  Gaussian 
noise. 

Let  ym(|  denote  a  sample  of  the  envelope  of  the  received  narrowband  proc¬ 
ess  on  the  m-th  alternative  signal  location  and  the  d-th  diversity  branch  (after 
narrowband  filtering  to  the  expected  signal  bandwidth),  where  1  <  m  £  M  and 
1  <_  d  <  D.  It  is  presumed  that  all  MD  random  variables  |ynid|  are  independ¬ 
ent  of  each  other.  Then  for  the  observation  vector 

x  =  JVA1 ...  y1D ...  y mi  •••  vmd]  • 

the  a  posteriori  probability  of  alternative  k  (after  observation  oi  z)  is 

Prob  |l4  p  (J  |  k) 

Prob  | k  | zf  =  * - ,  0  <  k  £  M  ,  (17) 

P(x) 
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where  alternative  k  -  0  denotes  no  signal,  Prob  |k(  is  the  a  priori  proba¬ 
bility  of  alternative  k,  and  p(z|k)  and  p(z)  are  the  conditional  and  uncon¬ 
ditional  probability  density  functions,  respectively,  of  z.  The  denominator  of 
(17)  can  be  expressed  as 

M 

P(Z)  =  ^  Prob  |  k }  p(z  |  k)  (18) 

k  =  0 


if  desired. 


Under  the  conditions  described  above,  one  has 


M  D 


p(xio>=  jj  n 

m  =  l  d  =  1 


ymd 


exp  - 


n 


__md 
2  a" 


,  y  .  >  o  , 
md 


(19) 


ykd 


2  2 

<r  +■  ir 

s  n 


exp  - 


2<% + 


(20) 


2 

where  <r  is  the  common  noise  power  level  in  all  MD  branches.  And,  for 
k  >  1,  n 

M  D  y  ,  /  ymd\  ^ . 

^■nn-r-y  n 

m=  1  d  =  1  %  \  2  n/  d=l 

m /k 

2 

where  <rs  is  the  common  signal  power  level  in  the  D  occupied  branches.  A 
more  convenient  form  of  (20)  that  is  available  is  (by  use  of  (19)) 


p(z|k)  =  P(z  1 0) 


'n 


2  2 

L<r  +  <r 

s  n< 


exp 


D 


2  2 


*8  +  %  d=  1 


ykd 


,  k  >  1  .  (21) 


Now  it  is  presumed  that  the  a  priori  probabilities  in  (17)  satisfy  the  follow¬ 
ing  rules: 


1-Q 

Prob  lol  =  Q ;  Prob  \k\  = - ,  l<k<M;  (22) 

'  ’  M  “  ” 


that  is,  the  a  priori  probabilities  of  all  the  signal  alternatives  are  equal.  Then 
selection  of  the  maximum  a  posteriori  probability  in  (17)  is  tantamount  to  pick¬ 
ing  the  maximum: 
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max 

0  <  k  <  M 


Prob  ]  k  |  z  |  corresponds  to 


But  this  is  equivalent  to 


D 

max  ^  2  > 

1  < k  <  m  2-j  Ykd  <  A  ’ 
d=  1 


(24) 


where  threshold  A  absorbs  the  a  priori  probabilities,  the  signal  and  noise 
levels,  and  D  and  M.  Satisfaction  of  the  upper  inequality  in  (24)  leads  to  a 
decision  of  signal  present  on  the  particular  alternative  that  dominates  the  left 
side  of  (24);  satisfaction  of  the  lower  inequality  in  (24)  leads  to  the  statement 
that  no  signal  is  present  in  any  of  the  M  alternative  locations.  *  Physically, 
(24)  indicates  that  the  squared  envelopes  on  all  D  diversity  branches  for  a 
particular  alternative  should  be  summed,  and  the  largest  sum  (out  of  the  M 
alternatives)  should  be  compared  with  a  threshold.  These  sums  are  the  deci¬ 
sion  variables  mentioned  in  (1)  through  (5).  For  a  specified  false-alarm  prob¬ 
ability  Pjta,  only  <r^  need  be  known  in  order  to  assign  a  value  to  A;<j-§  need 
not  be  known  (but,  of  course,  the  detection  probabilities  do  depend  on  <ri?). 

It  should  be  noted  that  if  the  cost  of  mistaking  signal  location  k  for  signal 
location  j  is  independent  of  k  and  j,  and  if  the  cost  of  mistaking  signal  loca¬ 
tion  k  for  no  signal  is  independent  of  k,  the  minimum  average  cost  results 
when  the  identical  test  (24)  is  utilized.  Of  course,  threshold  A  then  involves 
these  costs  also. 


A  second  signal  and  noise  model  that  leads  to  the  same  results  as  those 
above  is  that  of  envelopes  of  narrowband  deterministic  (except  for  phase)  sig¬ 
nals  subject  to  slow  Rayleigh  fading  and  additive  narrowband  Gaussian  noise. 
For  example,  |ymcj}  could  be  samples  of  the  envelopes  of  the  outputs  of  nar¬ 
rowband  matched  filters  that  are  nonoverlapping  in  time  or  nonoverlapping  in 

*This  derivation  parallels  that  in  reference  1,  section  IV,  very  closely. 
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frequency.  The  signal  strengths  in  the  D  branches  are  assumed  to  fade  inde¬ 
pendently,  and  the  signal  strength  changes  only  slightly  in  a  time  interval  equal 
to  a  narrowband-filter  impulse-response  duration.  The  quantity  <r|  is  the 
matched-filter  output  signal  power  per  branch  (averaged  over  the  Rayleigh 
fading  statistics).  The  relevant  statistics  are  identical  to  those  given  in  (19) 
and  (20)  [see,  for  example,  reference  4,  section  7.  5]. 

A  third  signal  model  is  afforded  by  samples  of  a  zero-mean  Gaussian  proc¬ 
ess  (rather  than  by  envelopes).  If  the  number  of  samples  on  each  alternative  is 
even,  such  as  2D,  then  it  may  be  shown  that  the  maximum  a  posteriori  proba¬ 
bility  processor  takes  on  exactly  the  form  of  the  test  in  (24),  where  each 
is  equal  to  a  sum  of  squares  of  two  zero-mean  Gaussian  variates. 

Thus  test  (24)  is  appropriate  for  three  different  signal  and  noise  models; 
its  performance  is  determined  in  the  next  section. 


PERFORMANCE  RESULTS 


The  false-alarm  and  correct-decision  probabilities  are  derived  in  appendix 


B  and  are  given  by 


PFA  1 


1  -  exp  (-A)  AkA! 

*-  k  =  0  J 


M 


(25) 


and 


CD 


I 


CO  X 

dx  — - 


D-  1 


exp  f  - 


(D-l)  !(1+R) 


D 


1  +  R 


D-l 

1  -  exp  (-x) 

S  xk/k! 

k  =  0 


M-l 

.(26) 


where  the  signal-to-noise  ratio  per  diversity  branch  is 

2 

ffs 

R  = -  . 


(27) 


n 


[The  result  in  (26)  is  a  generalization  of  equation  (19)  in  reference  1  in  that  the 
lower  limit  of  integration  is  not  zero.] 
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The  parameter  R  in  (27)  can  be  interpreted  differently  for  the  second  sig¬ 
nal  model  of  slow  Rayleigh  fading  signals.  For  matched  filtering  of  each  of  the 
D  diversity  branches,  it  is  shown  in  appendix  C  that 


1-1  E^/Nq 

N0  D 


(28) 


where  Ej  is  the  average  received  signal  energy  on  one  branch,  E-p  is  the 
average  total  received  signal  energy  on  all  D  branches,  and  N0  is  the  single¬ 
sided  noise  power  density  level.  The  curves  to  follow  will  employ  the  signal- 
energy-to-noise-density- level  ratio  (ENR)  parameters  in  (28). 


In  figures  1  through  35,*  the  probability  of  correct-decision  Pcd  is 
plotted  versus  E]/N0  in  dB,t  with  diversity  D  as  a  parameter  (solid  curves), 
for  M  =  1,  2,  4,  1G,  G4,  25G,  1024  and  PFA  =  10"n,  n  -  2,  3,  4,  G,  8.  The 
dashed  curves  on  each  figure  connect  points  of  equal  total  signal-energy-to-noise- 
density-level  ratio  Ep/No-  The  first  point  to  observe  from  these  detection 
characteristics  is  that,  for  a  given  available  Ep/N0,  there  is  an  optimum 
order  of  diversity  to  attain  the  maximum  value  of  Pqq  (for  a  specified  M 
and  PfA),  as  exemplified  by  the  peaks  of  the  dashed  curves.  However,  the 
maximum  PFD  is  not  always  realized  when  the  ENR  per  branch  is  5  dB,  2 
but,  rather,  it  can  be  significantly  larger  than  this  value.  For  example,  in 
figure  5  for  PFA  -  10-8,  M  =  1,  the  optimum  order  of  diversity  for  ET/N0 - 
13  dB  is  D  =  1;  therefore  the  optimum  ENR  per  branch  is  13  dB,  realizing  a 
value  of  PFF=  0.42.  And  in  figure  35  for  PpA  =  10“®,  M  =  1024,  the  same 
behavior  occurs  for  Et/N0  =  14  dB,  yielding  P£D  =  0.38. 

The  feature  of  large  ENRs  per  branch  being  optimum  is  characteristic  of 
the  smaller  values  of  Pcd  (that  is,  P^d  <  0.5).  For  larger  values  of  ET/N0, 
where  values  of  PFD  near  unity  are  attainable,  the  optimum  ENR  per  branch 
more  nearly  approaches  5  dB,  and  the  optimum  order  of  diversity  increases. 
However,  in  figure  35  for  PFA  “  10“®,  even  as  large  a  value  of  PFF  as  0.  9999 
still  requires  an  optimum  ENR  per  branch  of  E]/N0  =  7  dB,  not  5  dB.  Only 
when  the  two  "transition"  probabilities  —  (1  -  P^d)  PFA  (8)  (9)» 

respectively  —  are  approximately  equal  does  the  optimum  ENR  per  branch 
approach  5  dB.  But  for  large  M  (such  as  1024)  the  optimum  ENR  per  branch 
is  still  larger  than  5  dll,  even  for  1  -  P^q  ss  PfA- 


*All  figures  are  grouped  together  beginning  on  page  13. 
^That  is,  10  log10  (Ej/N0). 


9 


In  figures  36  through  38,  the  minimum  required  values  of  E'p/NQ  in  dB 
are  plotted  versus  M,  with  Pp^  as  a  parameter,  for  PcD  =  0.  5,  0.9,  and 
0.99.  The  results  are  extracted  from  figures  1  through  35;  results  for  other 
values  of  Pcd>  such  as  0.9999,  can  also  be  obtained.  The  slight  irregular¬ 
ities  in  the  curves  are  caused  by  an  inability  to  determine  the  dB  levels  from 
figures  1  through  35  more  accurately  than  0. 1  dB.  The  integer  numbers  under 
each  curve  denote  the  approximate  optimum  order  of  diversity  at  which  the 
minimum  E-^/Nq  is  realized. 

One  of  the  most  striking  features  of  figures  36  through  38  is  the  slow  rate 
of  increase  with  M  of  the  minimum  Et/Nq  required.  (It  must  be  remembered 
that  M  is  the  number  of  alternative  signal  locations,  not  the  number  of  signal 
observations;  this  latter  quantity  is  represented  by  D,  the  order  of  diversity. ) 
For  example,  at  the  high-quality  performance  level  of  PF^  =  10"8,  Pcd  = 
a  1  dB  increase  suffices  over  the  range  of  M  from  1  to  1024.  At  the  more 
moderate  performance  level  of  Ppy\  -  10-4,  Pcd  =  0.  5,  a  2.  5  dB  increase  in 
Et/N0  is  required  over  the  same  range  of  M. 

The  optimum  order  of  diversity  depends  very  strongly  on  the  desired  value 
of  PCd  but  is  somewhat  less  dependent  on  the  specified  PpA-  Figures  36 
through  38  show  that  the  optimum  diversity  is  1  to  2  for  Pcd  =  5»  5  to  10  for 

Pcd  =  0.  9,  and  11  to  19  for  PCD  =  0. 99.  If  the  optimum  order  of  diversity  is 
not  attainable,  due  perhaps  to  limited  available  bandwidth,  figures  1  through  35 
indicate  the  additional  amount  of  signal  energy  necessary  to  realize  the  speci¬ 
fied  level  of  performance. 


A  sample  program  for  the  generation  of  the  detection  characteristics  in 
figures  1  through  35  is  furnished  in  appendix  D.  Extension  to  other  ranges  of 
ppA,  D,  or  M  is  possible  by  modification  of  this  program. 


DISCUSSION 


For  nonfluctuating  signals  received  under  phase-incoherent  conditions  and 
with  D  =  1,  *  the  probability  of  correct  decision  is  given  [reference  5,  (78), 
with  crosscorrelation  coefficient  \=0]  by 


PCD 


dx  x  exp 


(29) 


*No  diversity  is  needed  for  nonfluctuating  signals. 
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where 


f2E 


(30) 


and  E  is  the  received  signal  energy.  The  false-alarm  probability  is  [reference 
5,  (77)] 


FA 


=  1 


exp 


(-  A2/2)] 


M 


(31) 


(These  results  are  identical  in  form  to  (14)  and  (11),  respectively;  only  the  den¬ 
sities  Px(- )  and  p0(-)  have  changed. )  Now,  (31)  could  be  solved  explicitly 
for  A  in  terms  of  PFA  and  M,  and  then  (29)  could  be  numerically  integrated 
to  evaluate  Pcd-  Alternatively,  the  bracketed  term  in  (29)  can  be  expanded  in 
a  binomial  series,  yielding 

pcd  -  r,  (~)X)  i-1*"-'  C)  e*p  (~)  <32> 

where  Q(-,-)  is  the  Q  function  [reference  4,  (4-55)].  However,  the  alter¬ 
nating  sequence  of  large  numbers  in  (32)  will  suffer  loss  of  significance,  and 
(32)  may  be  useless  for  large  M.  As  a  final  alternative,  the  upper  and  lower 
bounds  of  (A- 10)  could  be  used  for  evaluation  of  Pcd* 

A  comparison  of  the  performance  of  the  nonfluctuating  signal  results,  (29) 
through  (32),  with  the  present  results  has  not  been  pursued  any  further  except 
for  M  =  1.  In  this  case,  (31)  and  (32)  reduce  to 

PpA  =  exp  (-A2/2)  (33) 

and 

PCD  =  Q«i,A)  =  Q^, 


If  these  results  are  superposed  on  figures  1  through  5,  where  the  abscissa  is 
interpreted  as  10  logio  (E/N0),  it  is  found  that  (33)  cuts  across  the  D  =  1 
curve  at  approximately  Pcd  =  25  for  PpA  =  10-2  and  cuts  across  at  ap¬ 
proximately  Pcd=0.  36  for  PpA  ~  10”®.  That  is,  the  performance  for 


11 


TR  4793 


fluctuating  signals  is  better  than  for  nonfluctuating  signals  for  small  correct- 
decision  probabilities  (<  0.4);  the  occasionally  large  signal  amplitudes  occur¬ 
ring  in  fading  are  actually  helpful  when  the  noise  level  is  high.  Of  (fourse,  for 
larger  ENR,  where  PCD  approaches  1,  the  nonfluctuating  signal  case  yields 
better  performance  by  several  dB.  For  example,  Pp^  =  10-8,  PCD  -  0.99 
requires  an  ENR  of  15.4  dB  for  no  fading  versus  19. 1  dB  for  fading,  a  differ¬ 
ence  of  3.  7  dB.  The  same  difference  also  prevails  approximately  at  PpA=  10”2- 
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Figure  4.  Detection  Characteristics  for  Pj.-A-  10"6,  M  -  1 
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Figure  9.  Detection  Characteristics  for  P  ^  -  10  M  -  2 
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Figure  10.  Detection  Characteristics  for  PFA 
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Figure  13.  Detection  Characteristics  for  -  10-4,  M  =  4 
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Figure  14.  Detection  Characteristics  i'or  =  10  ,  M  =  4 


26 


TR  17!) 


E./N  (dB) 

I  o 


Figure  15.  Detection  Characteristics  for  P ^  -  10“y,  M  -  4 
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Figure  16.  Detection  Characteristics  for  PF» 
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Figure  17.  Detection  Characteristics  for  Pj  ^  =  10“y,  M  =  16 
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Figure  19.  Detection  Characteristics  for  PpA 
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Figure  21.  Detection  Characteristics  for  =  10  “2,  M  =  64 
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Figure  22.  Detection  Characteristics  for  PpA  =  10 
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Figure  23.  Detection  Characteristics  for  PpA  =  10”^,  M  =  G4 
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Figure  25.  Detection  Characteristics  for  P„.  =  10-8  M  =  64 
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Figure  26.  Detection  Characteristics  for  P^  =  10-2,  M  =  256 
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Figure  29,  Detection  Characteristics  for  P  .  =  10“6,  M  =  256 
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Figure  31.  Detection  Characteristics  for  P  =  M  =  .1024 


43 


Bim 

Ifj 

m 

MM 

pKl 

1  w!wiwi%wjk 

WffflMm 

MMM 

mmm 

f/SA 

YSA\ 


e  32.  Dete 


24 


0.9999 


TR  4793 


o 

u 

a. 


0.9999 

0.9995 

0.999 

0.998 

0.995 

0.99 

0.98 

0.95 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0.05 

0.02 

0.01 

0.CO5 

0.002 

0.001 

0.0005 

0.0001 


E,/N  (dB) 
I  o 


Figure  34.  Detection  Characteristics  for  =  10“6,  M  -  1024 
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Appendix  A 

ERROR  ANALYSIS  AND  INTERRELATIONSHIPS 
ERROR  ANA  LYSIS 


The  probability  of  correct  decision,  which  is  given  in  (14)  of  the  main  text 


as 


pcd‘[  ^PiMpy*)]"-1. 


(A-l) 


will  be  approximated  by 


A 

^  r  Ml  f°° 

PCD  =  JA  ^  Pl(x>  [p0(x)3  +j  dxp1(x). 


(A-2) 


The  error  is  then 


/•OO 

0  <  Error  =  PCD  -  PCD  -J  dx  p1(x)  |l  -  [P0(x)]M_1( 


<|l-[P0<A>]M'1}|l-P1<A)} 


(A-3) 


Now  we  must  choose  A  large  enough  that  Error  is  sufficiently  small.  One 
way  of  accomplishing  this  is  to  solve  for  A0  and  Ai  in 


1  -  [P0(Ao>]M  1  =  Error 

1  -  P  (Aj)  =  Error 


(A-4) 


and  then  choose  A  =  minimum  (A0,  A}).  Only  the  first  integral  in  (A-2)  need 
be  numerically  evaluated  inasmuch  as  the  second  integral  is  given  by  1  -  Pj(A). 


INTERRE  LA  TIONSHIPS 


From  (11)  and  (12),  we  have 


P  <  1  -  P 
MD  “  *A 


(A-5) 
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This  is  a  very  weak  inequality  in  most  practical  cases.  However,  in  the  limit 
of  small  signal-to-noise  ratios,  1  -  PMD  approaches  PpA* 

From  (13),  we  have 

M-l  M-l(  M)  M-l 

PID-,M'1){dXPoW[P<>Wl  =7|‘‘[V]  ’fa>  <A'(,> 

that  is,  the  incorrect-decision  probability  is  always  less  than  the  false-alarm 
probability. 

Upper  and  lower  bounds  on  PCD  are  also  attainable:  from  (5)  we  have 

^  PCD— Pro13  {X1  >  A  t S+N }  =  1  “  pi(A)  »  (A_7) 

whereas  from  (14)  we  have 

M-l  M-l 

PCD-  tP0<A)]  JA  dxPx(x)  =  [p0(A)l  C1  “  Pi^  •  <A_8) 

Therefore  it  holds  that 

M-l 

[p0(A)]  [1  -  px(A)]  <  pCD  <  1  -  pj(A)  •  (A-9) 

By  employing  (11),  this  can  be  expressed  as 

M-l 

t1  -  pfa1  M  t1  -  P1(A>] 

which  is  a  very  useful  and  tight  bound  for 
requires  no  numerical  integrations  at  all. 


-  PCD-  1  "  P1(A)  ’  (A"10) 

p  «  1.  In  fact,  the  use  of  (A- 10) 
FA 
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Appendix  B 

DERIVATION  OF  FALSE- ALARM  AND  CORRECT-DECISION 

PROBABILITIES 


From  the  section  on  definitions  of  error  probabilities  and  equation  (Z4),  we 
see  that  the  M  decision  variables  can  be  expressed  as  (proportional  to) 


1  <  k  <  M  . 


(B-l) 


When  the  probability  density  function  for  noise-alone  given  in  (19)  and  (20)  is 
utilized,  the  probability  density  function  pQ(* )  in  (11)  through  (14)  is  given  by 

x^1  exp  (-x) 

P  (x)  - -  ,  x  >  0  .  (B-2) 

(D-l)  ! 


Similarly,  the  probability  density  function  p,  (• )  is  available  from  (20)  as 


p  (x)  =- - - - -  ,  x  >  0  ,  (B-3) 

(D-l)  I  (1+R)D 


2  2 

where  R-ffg/o^  is  the  signal-to-noise  ratio  per  diversity  branch.  Then  the 
cumulative  distribution  functions  are 

D-l 

PQ(x)  =  1  -  exp  (-x)  ^  ^  xVk  !  ,  x  >  0  (B-4) 

k  -  0 


Pj(x)  =  1  -  exp 


(  X  1  /x  \ 

E  TT  “  *x>0' 


(B-5) 


The  false-alarm  and  assorted  decision  probabilities  follow  upon  substitution  of 
(B-2)  through  (B-5)  in  (11)  through  (14). 
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Appendix  C 

RELATIONSHIP  OF  SIGNAL^  TO- NOISE  RATIO  TO  SIGNAL-ENERGY  - 
TO-NOISE-DENSITY- LEVEL  RATIO 

For  transmitted  signal  R(t)  cos  [u0t  +  0(tjj  ,  the  received  waveform  on 
one  diversity  branch  in  slow  fading  is  given  by 

w(t)  =  A  R(t)  cos  [wQt  +  0  (t)  +  4>]  +  n  (t)  ,  (C-l) 

where  A  and  $  are  the  instantaneous  amplitude  and  phase  of  the  fading.  The 
output  cf  a  synchronized  matched  filter  (prior  to  envelope  detection)  is 


y  ~  /dt  R(t)  cos  [wot  +  ©(t)]  w(t)  =  yg  +  yn  .  (C-2) 


where 


1  2 

y  ^  —  A  cos  ♦  Jdt  R  (t) 
2 


(C-3) 


and 


y  =  Jdt  R(t)  cos  fw  t  +  0(t)l  n(t)  .  (C-4) 

n  L  o  J 


Then,  for  uniform  signal  phase  <|>  and  white  noise  n(t),  filter  output  powers 
are 


1 

4 


A2|-[JdtR2(t)] 


2 


No  i 
2  2 


Jdt  R2(t)  . 


(C-5) 


Now,  the  instantaneous  received  signal  energy  per  diversity  branch  is,  from 

(C-l), 

E^-A2  Jdt  R2(t)  .  (C-G) 

Therefore  it  holds  that 

EX^A2  Jdt  R2(t)  .  (C-7) 
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(C-8) 


By  combining  (C-5)  and  (07), 


one  obtains 


R 


[See  also  reference  1,  (22).] 

Notice  that  nowhere  in  the  above  derivation  has  it  been  necessary  to  specify 
the  time-bandwidth  product  of  the  component  pulse;  therefore  (C-8)  is  applicable 
to  any  time-bandwidth  product,  provided  the  fading  is  slow  and  is  frequency-non- 
selective  relative  to  the  pulse  duration  and  bandwidth. 
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Appendix  D 

PROGRAM  FOR  GENERATION  OF  DETECTION  CHARACTERISTICS 


DOUBLE  PRECISION  F.V1.  V2 
INTEGER  0(1?) *  D1 
EXTERNAL  F 

DIMENSION  Z  <  200  t  *  G1 12  I  t F  AC  1 12 1  » PR 08( 23.12) .TN0RMI25) 

COMMON  01 

DATA  !V  1.2.3.4.6.8.12.16.24.32.48.64/ 

DAT  A  YNORM/-  3.71902  »-  3  .2 90 5 3* -  3 . 09 023 »- 2 . 8 7 9 1 6 .-2 . 57 58 3 1 -2 . 326 3S * 
4-2.053 75.-1  .64485.-1  .  28 156. -.B 4162.-. 5 2 44  0. -.253  35.0..  .25 335*. 6 24 4 
*..B  A162*1.28 155.1  ,6<*<*  8  5 . 2 . 0537  5 . 2 .326  35  *  2. 5  7  5  S3 . 2. 3  78 16 . 3. 09023  . 
43,29053*3. 71902/ 

DO  21  IPF=2.8.2 
PF=10. **( -IPF 1 
DO  22  IM  =2 . 4  .2 
M=T  M 
M 1 -M-  1 
OMd./M 
FACtl 1=0. 

FAC  (2) =0. 

FAC  «3) =LOG<  2 .  I 
FACC4 1=1 0G(6.I 
F  AC ( 5 ) =  LOG( 120.) 

FAC<  6>=L  0CI5040.) 

F  A  C I  7 ) =  LOG(  3.99168E7  ) 

FAC<  8 1 =L 0GI1  .307674 37 E12) 

FAC(9)-L0G(  2. 5852016  7E 22) 

FAC( 101 =L  0GI8.22  2  83865E3  3) 

FAC<11  I=10G(2.5B623242I*59.*LCGU0.) 

FAC(12)-L0G(1. 982608321*87. *LCG(10.) 

S  =  0. 

T=-l. 

00  1  N=1.100 
T=T» ( N-l . -0M1 *PF/N 
S  =  S*T 

1  IF  <ABS( T).LT  .1  ,E-9»ABS( S)  )  GO  TO  2 

2  G(ll=-lOG(SI 
DO  3  30=2.12 
01=D(TD)-1 

Vl  =2 . »  G( 10-1 )  ♦  1 5  • 

V  2  =  VI* 1 .00 

CALL  ROOK  F  *OBLE(  S  ) .  VI .  V2  I 
Ct ID) =  V  2 

3  PRINT  4*  DIIDI.OCIDI 

4  FORMAT!'  D=*iX3t'  THRE SHO LD = • . E 1 5 . 9 > 

00  5  30=1.12 

D1=D(ID)-1 
DO  6  IE  =  1 . 23 
ETNO  DB  =IE ♦ 3 

R=10.** I . 1»ETN008 I/D ( ID  I 

0=1. n  i.*r) 

Q=(  l.*R  1  •  t  3  3.  *1 .38*01  ID) -15.*  EX  PI-  .092*01)  ) 
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COKP-DI ID  I »L OG( 1 .♦H) »rflC  I  ID I 
u :S  IMPS  1 H » 51 ID) .U . I23> 

PRINT  7.  01  ID  I »  E  TNOD S . V 

7  rORMATII10»E10.2.£20.81 
V=SIMPS2  10) 

PRINT  8 » V 

8  F0RMATIE40.3) 

V  =  MINl V». 9995^9) 

PROBIIC«ID)=HAX<V.l.E-o> 

6  CONTINUE 

5  CONTINUE 

CALL  NODE  SC  I  2 ■ Q ) 

CAL  L  SU9J  EG ( l . 0.. YNOR  Mil )  *  2C. . Y NOR  HI  25) ) 
CALI  OB JC  T  C I Z >1150.1  ICO  .  »  285D.i2700. ) 
CALL  StTSMCIZ.30.1.) 

DO  9  J  -  2  >  1  9  »  2 

CALL  LINESCI2.C.FL0AI(JI.YN0RM<1)I 

9  CALL  l  INESGfZ.l.FLOATt  J|  .YN0RHI25H 
DO  10  J=2.24 

CALL  LINE5GIZ.0.  O..YNORMIJ)) 

1C  CALL  LINE5DIZ.1.2D.. YNCRM(J)  1 
CALL  SETSMGlZ.30.2.) 

CALL  LTNESCIZ.O.  C  .  .  YNCRH  11)1 
CALL  L  INESG  (  Z  *  1  *  0  .  *  Y  NOR  Ml  2  51  ) 

CALL  LINESCIZtl.ZD.i YN0RMI25  I) 

CALL  L  INE5GIZ.1.20,  rYNORMID) 

CALL  HNE5GIZ.lt  C..YNORMUIJ 
CALL  SETSMGIZ  .8*4  .  »4»  ) 

DO  11  10=1.12 
DDB=10.*l  OGIOID  I  ID  I  ) 

P=TINORMIPR09I1.ID)  .1  3Q) 

CALL  LINESGIZ.0.4.-DD3. P  I 
CALL  POINTCIZ.1.4.-OD6.P) 

DO  11  I£  =  2. 2  3 
P=TIN0RMIPR0B(IE.ID).*3Q) 

CALL  LINESGIZ.:.IE-»3-DD8.P) 

11  CALL  POINTGIZ. 1. IE.3-DDB.P) 

CALL  PAGECI Z.3.4 .1) 

22  CONTINUE 

21  CONTINUE 

CALL  EXITGIZI 
3C  STOP 

EL'NCTTCN  HI  XI 
H  =  EXP  l-X) 

IE  101. EC. D I  GO  TO  2 
T  =  H 

CO  1  K  =  1 *  D  1 
T  =  T  •  X /K 

1  H=H ♦ T 

2  H  =  EXP (D1»L  OG IXJ-0 • X »M1 «  LOG  1 1 H) -  COMP  ) 
RETURN 

E  ND 
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SUBROUTINE  R00TIF.D.V1.V2) 

D0U8L  E  PRECISION  F  » O*  V  1 «  V2  .  X  (  1 0  Oi  .  OA  (1  OQ  )  »B  »B  A 
X  C 1 )  -  V  I 
Xt  2  >  =  V2 
DA  (1  >  =  F<X ( 1 >  ) 

DA<  2)  =F  I  X  f  2  t  ) 

N1  :2 

DO  9  N=3.10a 

TFUBSIDA(N-1)-D1.LE  .1.0-  9*ABSI0>1  GO  TO  5 
IF (DAI  N-l >-DA(N-2l 1  7.6*7 

6  X(NI=X(N-1)»1.100 
GO  TG  B 

7  Xt  N  );  (X  IN  -  1 1  *  C  0-  DA  IN-  2)>  +X  C  N-  2)  •  ( OA  I  N-  1  )-  0  )  )  /  <  DA  I  N-l  I  -OAt  N-2  )  ) 
B=X(N)-X<  N-l I 

BAiABSt  B ) 

XINI =X ( N-l )*MIN<0A.X (N-l I*. IDO) «SIGN( 1 .DO.B  » 

XI N  >  =  MAX (X  tN  » .0 .001 

a  ni  =n 

R  DAI  N>  =f  t  X (N I  J 

S  PRINT  9  •  XINI I.DAINl  I » N 1 

9  FORMAT  12 2D3C.18. 110) 

Y2  =  X  INI ) 

return 

END 


FUNCTION  SIMPStF.Xl.X2.NNI 
DOUBLE  PRECISION  ACE.ACO 
N  =( (MAX  I ABSINN  I  .1  I +1 ) 221 »2 
H-  (X2-X1 1 2 N 
D-2.*H/3. 

S  =  .S«( FI  XI ) ♦ F ( X 2  )  ) 

A  CE -0  .DO 

IF  tN.FQ.2 I  GO  TO  3 
LTM=N-2 

DO  2  T=2. LI N.2 
T=FIX1+H*II 

2  ACE=ACE*T 

3  A  CO -0 . DO 
L  IM=N-1 

DO  1  1  =  1. LI M. 2 
T=FIX1*H*I) 

1  AC  0  =  AC  0*  T 

S  TMi1  S  =  !  j  ♦  2  •  *  SN  Gl  I  ACO  )  ♦SNGLIACE)  J«0 
RETURN 

ENTRT  SIMP  S  2  1 1 DUMM  Y I 

NAN*  2 

H=H* . 5 

t>=  D*  .  5 

ACE  =  A  CE  *A  CO 

GO  TO  3 

END 
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DOUBLE  PRECISION  FUNCTION  F<X) 
DOUBLE  PRECISION  X»T 
INTEGER  02 
COKMOM  D1 
F  :E  XP  I  -  X  I 

I F I  01 .EQ.O)  RETURN 
T=r 

DO  1  K=1.D1 
T :T • XX  K 
F  =  F+T 
RETURN 
ENO 
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